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1. Introduction 

The staircase model of Al. Zamolodchikov is a simple relativistic factorized scat- 
tering theory in 1+1 dimensions, which shows signs of a very non-trivial renormalisation 
group behaviour. It describes a single boson with mass m and two-particle scattering 
amplitude 

SW = tanh(i^-^)ta„h(^-^^), (1.1) 

where Oq is a real parameter. Assuming the existence of an underlying field theory, the 
model can be studied at all length-scales by placing it on a cylinder and varying the 
circumference R. In particular, the ground-state scaling function c{x) {x = log can 
be obtained by means of the Thermodynamic Bethe Ansatz (TBA) (The ground- 

state scaling function can be interpreted as an 'effective central charge' for the non scale- 
invariant theory, and is related to the vacuum Casimir energy by E[R) = —nc{x)/6R; 
for a unitary scale-invariant theory it is just equal to the central charge.) Plotting c(x) 
as a function of x shows a series of plateaux, connected by steps each time x is near an 
integer or half-integer multiple of Oq. This staircase-like pattern becomes more pronounced 
as 9o increases, the values taken by c(x) on the plateaux then running through the series 
Cp = 1 — 6/p{p+l), the central charges of the c<l minimal models. This suggests an RG 
flow from an ultraviolet fixed point (i? — > 0; x ^ — oo) with c = 1, to a trivial c = fixed 
point in the infrared, passing close by each minimal model Aip in turn. Varying results 
in a one-parameter family of such fiows; the larger ^o, the more closely the trajectory visits 
each minimal model in its journey from c = 1 to c = 0. 

From another point of view, A. B. Zamolodchikov 0], and Ludwig and Cardy showed 
some time ago that for large p the deformation of Aip by its ^13 operator leads in the 
infrared to Aip-i, so long as the coupling constant is positive. More recently, TBA analysis 
by Al. Zamolodchikov [^,0 , applicable for all values of p, has reinforced this picture. The 
perturbed model, commonly denoted A4A^\ can thus be associated with an RG fiow that 
hops between two neighbouring minimal models, and the staircase model can be seen as an 
approximation to this whole series of hopping fiows, the approximation becoming better 
as 6*0 increases. 

The minimal models are only the first of many infinite series of conformal field theories 
that can be constructed by means of the GKO coset construction § . For example, given 
any simple simply-laced Lie algebra Q = A, D or E, there is a series of VFc;-minimal 
models, rational conformal field theories described by the coset Q^^'^xQ^'-^ /Q^^~^^\ Since 
the minimal models are recovered by chosing ^ = Ai, it is natural to ask whether staircase 
models can be found for the other Wg series. The fact that Zamolodchikov's S-matrix 
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( p. ■ 1|) is an analytic continuation in the coupling constant of the sinh-Gordon S-matrix 
(already remarked in ref. [Q) provides a strong hint that the required generalisation is to 
be found in the analytic continuation of the real-coupling afRne Toda S-matrices, and this 
turns out to be the case ■ To complete the analogy with the minimal sequence, it can 



be argued that 'one-hop' flows similar to those associated with the models A4A^'^ should 
also exist for each Wg series. The generalisations of Zamoldochikov's Ai TBA analysis [H 



to arbitrary ^ [11,12] lend support to this view. 

However this is by no means the end of the story. In particular, there are the 
models x^(')/^(^+') For each flxed k, they form a series with central charges 

Q = r{h+l)kl{k+2h+l) / {h+k){h+l){h+k+l) , where r is the rank of Q and h the Coxeter 
number. In all cases, the operator analogous to ^13 is (/)id,id,adj , the three indices labelling 
particular representations of G^'^\ Q^^^ and ^('^+^) respectively. 

For Q = Ai, there are general perturbative results for large / matching those already 
described for the minimal models There is a surprise here, in that the perturbation 
of the model of the k^^ series does not flow to the /— 1^*^ model of this series, but rather 
to the l—k^^. As a result, there is no longer a single sequence of one-hop trajectories, but 
rather k disconnected sequences, interlaced along the k^^ series. Note that taking Q = Ai, 
k = 2 gives the A^=l superconformal discrete series, for which this result had already been 
found in ref. [jl5|. A TBA analysis for the general Ai case, for arbitrary values of /, has 
now been given by AI. Zamolodchikov ; the picture outlined above holds good modulo 
a small complication when I becomes smaller than k, which will be described later. The 
further generalisation to arbitrary Q can be found in ref. [l^ . 



Since for given k there are k different sequences of hopping flows, it is reasonable to 
hope for k different staircase models approximating them. Some nice recent work by Mar- 
tins []T6| , p!7[| has started this programme by proposing a generalisation of Zamolodchikov's 
original staircase model to one of the k = 2 sequences, flnding TBA systems indicative 
of RG trajectories which pass close by the subset of superconformal minimal models with 
non-zero Witten index (/ even), before flowing off to massive theories in the IR limit. He 
also showed the generalisation of this to other Q, again at k = 2 and I even. One interesting 
feature of his proposal is that each TBA system contains magnonic terms (pseudoenergies 
with no direct couplings to energy terms), perhaps indicating that the as-yet unknown 
scattering theories underlying the models have non-diagonal S-matrices. 

In this paper, we propose and start to analyse TBA systems for general k and Q. The 
proposal itself is given in the next section: for each k there are k different systems, labelled 
by an index s G Z^. In sections 3 and 4, we describe why these systems should mimic the 
desired hopping behaviour, a discussion that has many parallels to one given earlier for 
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the k = 1 series [10 . The conclusions of these two sections have been backed up by some 
numerical work, which is reported in section 5. The RG flows predicted have a number 
of surprising features, which can however be supported by alternative arguments. These 
points, along with some further speculations, are contained in the concluding section. 



2. The spiral staircase 

Take Q to be an arbitrary simply-laced Lie algebra, with Coxeter number h, and 
let a label nodes on the corresponding Dynkin diagram. In real-coupling affine To da field 
theory, the different types of particle are labelled by just such an index; let nia be the mass 
of the corresponding particle. These masses can also be characterised as the components 



of the Perron- Frobenius eigenvector of the Q Cartan matrix [1T8H20|| . To provide for some 
magnonic structure, an extra label i will be needed, lying in Z^. As already mentioned, 
the particular staircase out of the k possibilities is determined by a member of Z^, s say, 
which will stay fixed throughout the discussion. 

The TBA system will be given in terms of r x k pseudoenergies ■* (r is the rank 
of Q). The Zfc-condition on the index i amounts to identifying Sa"^ with Sa^^K Similarly 

(i) 

there are r x k energy terms i^a , given by the formula 

zy«(6») = m„(5*'°e=^-^ + 5^'"e^+^) , (2.1) 

where 9 is the rapidity, and x encodes the scale of the system via = \miR. The case 
s=0 will turn out to have a massive infrared limit, allowing mi to be interpreted as the 
mass of a particle, or multiplet of particles; the values of the remaining nia are contained 
in the (dimensionless) ratios fha = ma /mi. Note that the energy term is only nonzero 
if i is equal to or s: it is in this way that the value of s enters the game. Defining 
La\o) = log(l + e~'^'°-^^^^), and introducing a parameter 6*0 G R"*", the proposed TBA 
system is: 

+ * - ^'^^ * ^^'V - ^o) - ^aB * lI:+'\0 + 6o)] = i^i'\e) , (2.2) 

b=i 

where all the i-type indices are to be taken modulo /c, and * denotes the rapidity convo- 
lution: (P * L{e) = J^^d9'(j){9 - 9')L{6'). The kernel funct ions (pab and ipab are defined 
in terms of the minimal parts of the corresponding affine Toda S-matrix elements 5"^'", 
written in the form 

ST{0)= n - 1) (^) + 1) (^) ; ix){9) + 
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where Aab is a set of integers (possibly with repetitions), and the related functions 



Then 



xeAab 



.6(ff) = -i^logS'"'»(») ; MO) 



More explanation can be found in refs. [|T2i|T0| ; while this notation is consistent with that 



of these two references, it differs from that of ref. [T^]. A complete list of the functions 
5'^^" can be found in ref. fl^, and general formulae in ref. ||21|| . 

The ground-state scaling function is now given in the standard way in terms of the 
solutions to the TBA system (|2.2|) , and the energy terms ( |2.1|) : 



c{x, 0o] 
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TT 



r fc — 1 

.EE 

a=l i=0 



rf^m„e^z.«(^)LW(^). 



(2.3) 



The next section will discuss the solutions to ( |2.2D , explaining how the scaling function 
c{x, 9o) can be expected to behave as a function of x. Before giving these details, it is 
worth making a few general comments on the system defined by ( ^.11 ) and ( |2.2D . First, 
notice that just as in refs. |]T6| , [T7[] , there are many 'magnonic' pseudoenergies with zero 
energy term: this may be a sign that for these models too, the underlying S-matrix is 
non-diagonal. More striking is the form that the shifted terms in ( p.2|) take: while the 
term shifted by —9o couples with L 



that shifted by +6*0 couples with L 



Thus 



these factors are not a simple continuation of the affine Toda Z-factors, as was the case for 
previous staircase TBA systems [plP|, [T0| , p!6t|T7| . Despite the apparent asymmetry that this 
implies, the solutions of ( p.2|) are symmetric under 9 —>■ —9, so long as this transformation 



is accompanied by an exchange of the pseudoenergies ^\ The possibility to 

implement the parity transformation in this way is a reflection of the symmetries of the 
affine Dynkin diagram on which the labels i live, and is reminiscent of the discussion 
for systems based on non-afhne graphs given in ref. p2|. It also allows the 'left-moving' 
and 'right- moving' energy terms in (|2.1| ) to be coupled to different pseudoenergies for st^O, 
a hint that in these cases the system has a massless infrared limit. More detailed analysis 
in a later section will confirm this expectation. 



The most important feature of (|2.2|) is the non-local nature of the interactions between 
the different pseudoenergies. The functions ^(6*) and ^^{9) are exponentially damped away 
from a region of order one about ^=0, so their convolutions 4) * L{9) and i(j * L{9) are, up 
to exponentially small corrections, dominated by the values of L{9') in a region of order 
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one about 6' =6. (While functions L(6) couldhe found for which this does not hold, suffice 
it to say that these do not seem to arise amongst the solutions to (pT^).) Hence the three 



convolution terms in (2^), which involve ea \ and e^*^^'*, pick up most of their values 

near 6*, 9—9q, and 9+9q respectively. To visualise this, it is helpful consider the index i on 
the same footing as 6*, albeit only taking discrete values. Thus we write ea\9) = ea{i,9), 
the pair (z, 9) being valued in Z^xR. This space has the form of a cylinder, on which ( |2.2| ) 
induces interactions between the neighbourhoods of (z,6'), (z— 1,6*— 6*0) and 9+9q). 
This explains the sobriquet 'spiral staircase': the pseudoenergies couple together in a 
spiral pattern round the cylinder, a feature that turns out to be important in reproducing 
the expected hops in the RG flow. 



3. The double helix 



We now turn to the solutions to and to progress we will have to make some 



more assumptions about their general form. These parallel the assumptions made in the 
analysis of more usual TEA systems, and we will not pretend to give any rigorous proofs. 
However, we have also made a number of numerical checks, which will be commented on 
later. 

Comparing ( p.2|) and (|2.1|) , there are two regions where Saih 9) and La{i, 9) are imme- 
diately known (as always, up to exponentially small corrections - we will tend to assume 
such a phrase to apply globally from now on). For i = and 9 <^ x, the energy term z^a(0, 9) 
becomes very large and dominates (|2.2|), so that £a(0,6'^x) ~ mae^~^ . Correspondingly, 
La(0, 9<^x) suffers a double-exponential decay and is soon negligably small. In the region 
i = s, 9 ^ —X, similar considerations show that ea{s,9^ — x) grows exponentially, while 
La(s, 9^ — x) quickly decays towards zero. 

Near to (0,x), all the terms in (|2.2|) come into near-equal competition, and it is no 
longer a good approximation to ignore the convolutions. However, beyond this transitional 
region the form of the equation simplifies once again, as for 9 ^ x the energy term z^a(0, 9) 
becomes exponentially small and can be dropped. On the basis of previously-studied TEA 
systems, we expect La{0,9) to have a kink near 9 = x, interpolating between zero for 
9 <^ X, and some other constant for 9 ^ x. Just as is the case for other TEA systems, 
the precise form of this kink is hard to find due to the increased complexity of (|2.2| ) in the 
transitional region, but is not required for an asymptotic evaluation of the ground-state 
scaling function. 

At this point the non-locality of ( p.2|) comes into play. The presence of a kink near 
(0, x) has an effect, via the shifted convolution terms, on the equations obeyed near 
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(/c— 1, a;— ^o) and (1, x+^o)- This is similar to the propagation of kinks in the k=l staircase 
model described in [1^, though with one important difference: here, the influence is on 



pseudoenergies with different values of i from that of the 'seed' kink near (0, x), and hence 
with different energy terms derived from ( |2.1| ). In particular, while La{0,9) was forced to 
be zero for ^ ^ x by the dominance of Ua{0,9) in this region, this is generally not the 
case for La{k—1,9). Thus we expect to find kinks generated in both directions from the 
initial kink. These in turn cause there to be kinks at (/c— 2, x—29q) and (2, x+29q), and so 
on, the set of secondary kinks spiralling round the cylinder in both directions from (0,x). 
Making the transformation i s—i, 9 —9 reveals another set of kinks spiraling away 
from (s, —x) , and interleaving with the first spiral to form the pattern of a double helix. 
This is illustrated in figure 1, where the kinks should be imagined to be strung out like 
beads along the two spirals. The two regions where the energy terms are dominant are 
depicted as double lines; the two seed kinks are at their ends. That they do not overlap 
means that the figure shows a situation where x is negative. 

The asymptotic directions for which c(x, ^o) can be evaluated correspond to the kinks 
from one spiral becoming far away from those of the other, so that the values of and 
La have time to settle down to approximately constant values in the inter-kink regions 
(that they do become nearly constant in these regions is in fact the key assumption of the 
discussion). This breaks down whenever {m,x+n9o) ~ (s, — x) for some m. The (strict!) 
equality implied on the Z^- valued first coordinate requires m = nk+s for some integer n, 
and so interleaving will fail each time x/^o ~ —{nk+s)/2 , corresponding to a crossover in 
the critical behaviour. Attention will now be restricted to situations far from these points, 
and the scaling function c(x,^) treated in the limit 9o ^ oo with x/9o remaining fixed at 
a value away from crossover. 

While secondary kinks are generated in both directions from each seed kink, they 
cannot propagate indefinitely. The fact that is cyclic means that at some stage energy 
terms will be re-encountered in domains where they are dominant, and this truncates the 
chain of kinks. For the remainder of this section, we assume that x < 0, as in figure 1; it 
will turn out that most (though not always all) of the hops are found in this region. For 
definiteness, take x to lie between crossovers in the region defined by 

-{nk + s)9o < 2x < -{{n-l)k + s)9o ■ (3.1) 

Consider the spiral generated by the kink at (0, x) (the behaviour of the other spiral 
will then follow by parity). After k—1 steps in the —^-direction, there is a kink near 
(1, X— (/c— l)^o)- The energy term Ua is zero in this region, and so there is no reason for 
this kink to be suppressed. However, one more step arrives at {0,x—k9o). This is well 
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inside the region i = 0^ 6 <^ x where, as aheady described, the part ma5'^'^e^~^ of the 
energy term dominates the TBA equation and forces La{0,9) to be vanishingly smaU. 
There is no kink here - La{0,9) has been zero since the seed kink at (0, x), and remains 
so as 9 decreases through x—k9o. Thus in the —^-direction, the chain of kinks terminates 
after k—1 steps. In the +9 direction, the truncation is instead effected by ma5^'^e^~^^ , the 
part of the energy term which couples when i = s. If ^ ^ —x this term plays no role in 
( p.2|) , and for x in the region (|3.1J ) the spiral makes n complete circuits of the cylinder, 
and then continues with a further s—1 kinks, before the energy term becomes important. 
So, there are nk+s—1 kinks in the +^-direction from the seed kink at (0, x), truncation 
being caused by the suppression of the putative kink at (s, x + {nk+s)9o). 

The full picture for both spirals can be seen by referring once again to figure 1: the 
spirals terminate when they encounter the double lines representing regions in which the 
energy term is completely dominant. Since they each make a single turn in one direction, 
and two complete turns in the other, the figure should correspond to taking n=2 in ( |3.1|) : 
this is easily verified by writing down the inequalities implied by the ordering of the marked 
points along the 6'-axis. 

We can now unwind each spiral from the cylinder, and represent the set of interacting 
kinks graphically: 

1 k-l k k+1 k+nk+s-1 

O O O O O (3.2) 

Each node is a kink, and nodes are linked if they interact via the shifted convolution terms 
in ( |2.2|) . Empty nodes ( Q ) represent kink regions where the energy term plays no part 
in the leading behaviour of (pT2|), while the filled node ( (3) ) corresponds to the seed kink, 
in the neighbourhood of which the energy term ma5*'*^e^~^ cannot be ignored. This is also 
the only region of those represented on the graph which contributes directly to the formula 
( p.3|) for c(a;, ^o)- An isomorphic (reflected) graph results for the other spiral. 

The next task is to find the values taken by the functions £a(z,^) and La{i,9) in 
the inter-kink regions. The necessary constraints follow in the usual way by pulling near- 
constant terms out of the convolutions in ( |2.2| ), leaving only overall integrals of (pab and 
tl^ab- These are[|g,g: 



- / d9<Pabi9) = 5ab - 2C-,i ; — / d9i;am = -C;,\ (3.3) 

where Cab is the (non-affine) Q Cartan matrix. The various inter-kink regions influ- 
ence each other through the non-local terms in (|2.2|), in much the same way as did the 
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kinks. To describe the situation, let the p kink in (|3.2|) be located near {ip, 9pj 



{pmodk,x+{p—k)9o), and let {ip, 9p ) and {ip, 9p) be points in the inter-kink regions im- 
mediately before, respectively after, this kink: 

9; = ^9o{2p-{n+2)k-s) ; 9^ = ^9o {2p - {n+l)k - s) . (3.4) 

Via ( |2.2[ ), the ea{ip, 9p) come into interaction, as do the ea{ip, 9p). However care is needed 
when p=k. While the seed kink at (ifc, 9k) = (0, x) owes its very existence to the delicate 



balance between the energy and convolution terms in (|2.2|) at that point, the energy term 
has become negligable by the time that {ik,9'^) is reached, and so has no effect on the 
equation for ea{ik^ ^k)- other hand, {ik, 9^) is in the region immediately to the left 

of the seed kink, where the energy term in ( p.2| ) is dominant. As a consequence ea{ik-, (^k) 
is forced to be effectively infinite, irrespective of the values taken by the other ea{ip, 9p). 

All this has the effect of changing the 'connectivity structure' that for the kink inter- 
actions was summarised in ( |3.2| ). It can be visualised on figure 1 by mentally shifting the 
spiral based at (0, x) slightly to the left, and then slightly to the right. Shifting leftwards, 
the spiral is cut by the double line representing a dominant energy term, while shifting 
rightwards disconnects the spiral from this line altogether, apart from the endpoint. Thus 
for the inter-kink regions containing the left-shifted points {ip, 9p), ( |3.2|) should be replaced 
by 

1 fc-1 fc+1 fc+nfc+s-l 

O O O O O (3.5) 

while to the right, the graph for the {ip, 6*+) is 

1 fc-1 fc fc+1 fc+nfc+s-l 

O O O O O O (3.6) 

For the second spiral making up the double helix, kinks are found at (Jp,9p) = {s — 
V' ~^p) the two neighbouring inter-kink regions are centred on (ip,9^), where 9p = 
—9~ = 97 and 9~ = —9j~ = 97 ..... . The above treatment goes through 

P {n-\-2)k-\-s—p P P (n+ljfc+s— p 

isomorphically when phrased in terms of these kinks, as indeed it must by the parity 
symmetry of the system. 

To express the algebraic content of ( p.5|) and ( p^) , let /"^ and Ip^ be their respective 
incidence matrices, and put (7^ = 25pq — Ip^. At {ik,9^), the value of has already 
been established; for the rest, extracting the near-constant values of ea{ip,9^) from the 
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convolutions in (|2.2| ) and using ( p.3|) yields, after some simple rearrangements, the following 
consistency conditions: 



r (n+i)fe+s-i n-^n 

' ^ab '-'PI 



where 



6=1 0=1 



These equations reproduce the constraints for the limiting values of the pseudoenergies 
in the Qi'^) y,gi'^^+'') /gi^'^+'^)^+'') TBA system given in providing the substitution 

exp(— £a(zp, 6'^)) = Y^i^^oo) is made. While this is suggestive, caution is needed, since 
( p77|) describes the values of the functions £„ at various points along the 6'-axis, rather than 
just at ±oo as for more usual systems. Nevertheless, asymptotically far from the crossovers 
the different kinks should become decoupled and behave like independent functions, their 
only interactions being those explicitly marked on ( |3.2|) . These reproduce exactly the 
interactions between different pseudoenergies found in the TBA systems of 0,0, so we 
certainly expect the evaluation of the function c{x^6q) in the asymptotic regime under 
discussion (given by (|3.1j)) to be the same as for the corresponding coset TBA system. But 
to confirm this, the leading asymptotic behaviour of c(x, ^o) away from crossover can also 



be calculated directly from the formula (|2.3| ) , remaining within the spiral staircase model. 



The necessary modifications to the usual arguments have already been described at some 
length in ref. and the slightly increased complexity of the models being considered 
here has no bearing on the calculations once the spirals have been unwound, and the 
connectivity structure ( |3.2| ) established. The one small point to watch is that since the 
massive kink is now located at the k^^ position, inside the chain, the 'integrations by parts' 
by means of which the other kinks are felt now proceed in both directions away from /c. 



rather than just to the right as for the k=l staircase. Referring to the earlier paper [ffD 
for further details, we will simply report the final result: 



„ r (n+l)fc+s-l 

D 



hin c{x, eo)=Cr. = -J2 E (/i) - ^ ifap)] , (3-9) 



where x/9q is held fixed at N^, —{nk+s)/2 <^ ^ — ((n— l)/c+s)/2 , while the limit is 
taken in order to stay inside the regime defined by ( |3.1j ). The functions appearing on the 
righthand side are Rogers dilogarithms, defined by C{z) = — \ dt _|_ ^ and 

their arguments are obtained by solving (|3.7|). To finish the calculation, a sum rule 
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is needed [Q: if {fap} i^^a^r, l<p<l) is the solution to ( p.7| ) when Cab is the Cartan 
matrix of Q and Cpq the Cartan matrix of the algebra Ar' , then 



TT 



a=l p=l 



rr'h' 
h + h' 



rh' - 



r{h+l)h' 
h + h' 



(3.10) 



where h is the Coxeter number of Q and h'=r'+l the Coxeter number of A^'. The second 
version is useful because it shows the sum to be rh' minus the central charge of the 
gih') wzw model. For the case in hand, the sum over the jC{f^p) follows from ( |3.1C| ) 
with h' = {n+l)k+s, since is the Dynkin diagram of . For the {f~p} 

there are two disconnected parts to (|3.5|) , the Dynkin diagrams of Ak-i and Ank+s-i- 
Correspondingly the dilogarithm sum splits into two, summed by taking h' equal to k and 
nk+s respectively. The absent k^^ node does not contribute since C{f~j^) = C{Q) = for 
all a, and so 

k nk+s {n+l)k+s 



r{h+l] 



+ 



h+k ' h+nk+s h+{n+l)k+s 
the central charge of the ^C^) x^("'''+'*)/^(("'+^)''+'^^ coset model. 



(3.11) 



4. The infrared limit 

To complete the picture of the asymptotics of the ground-state scaling function 
c{x,9q), the condition that x be negative, imposed thoughout the last section, must be 
relaxed. Since x is related to the physical parameters of the model by ^miR = e^, this 
corresponds to examining the system in the infrared. 

The point to note is that for x positive the nature of the relationship between the 
spirals of kinks and the two regions of energy-term dominance, illustrated in figure 1 for x 
negative, may change. When x becomes positive, the two regions, marked by double lines 
on figure 1, start to overlap, as in one region 9 runs from — oo to x, in the other from —x 
to -|-oo. If s = 0, then the two regions actually collide and this has the effect of killing 
off the two seed kinks completely. Thus in this case the IR limit is simple: a theory with 
c = 0, with the final crossover to a massive phase occurring near a; = 0. Otherwise, kinks 
will continue to be found near (0,a;) and (s, —x) for arbitrarily large positive values of x. 
Thus the Li'\e) never become trivial, and a non-zero value for c(x, ^o) is expected even 
as X +00. This reinforces the assertion made in section 2 that for s ^ the IR limit of 
( p.2|) should be massless. Precisely which conformal field theory this massless limit should 
be brings one last surprise. In the previous discussion it was seen that for x<0 the spiral 
generated at (0, x) truncates after a single turn in the —9 direction, when it re-encounters 
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the i=0 energy term. This is why the left-most piece of (|3.5| ) has k~l nodes, and one of 
the elements of the coset whose central charge ( p.9| ) reproduces has level k. But if x is 
positive, this can change. In the —9 direction, the chain from (0, x) might encounter the 
energy term at i=s before returning to i=0. The onset of this phenomenon is signalled 
by a crossover around x = ^{k—s)9o, beyond which point the chain of kinks anchored at 
(x, 0) is terminated in both directions by the i=s energy term. It is straightforwardly seen 
that this is the last such crossover expected on the basis of changes to the overall kink 
structure, and that after this point the situation stabilises, the form of the pseudoenergies 
remaining unchanged apart from simple translations as x ^ +oo. The final kink system, 
governing the IR limit, can be unwound and represented graphically just as before. The 
picture (|3.2D becomes: 



s+l fc-l k k+1 k+s-1 

cy-o — <s> — oo (4.1) 

where to ease comparison with the earlier graphs, the mod- A; values of the labels have 
been preserved. This graph has k—s—1 massless nodes to the left, s—1 to the right. The 



remaining calculations now go through unchanged, (|3.9|) becoming 



lim c(x,6'o) = c_i = r{h+l] 



"0- 

a:/eo=JV_i»(fc-s)/2 



k—s s k 

+ 



h+k—s h+s h+k 



(4.2) 



which is the central charge of the Q^'^~^^ xQ^^^ /Q^'^^ coset model. 

In this last equation, the fact that A^_i is positive means that x tends to -|-oo as the 
limit ^0 ^ oo is taken, but this does not spoil the validity of the asymptotic estimates 
being made. However, the real interest all along has been to trace the variation in c(x, ^o) 
in one particular model, for which is fixed. Thus a change in the point of view is needed 
to apply the results found so far. If the value of Oq is large enough, then as x varies, 
the pair (x, 6*0) will pass through a series of regions for which the asymptotic results ( ^?9| ) 
and (^3) are good approximations to the true values of c(x, 6*0). Hence this function will 
run through the values in turn, deviating significantly from these numbers only in the 
crossover regions. An order-of-magnitude estimate for the size of these regions is easily 
given: the approximations leading up to (^.9]) and ([4.2|) were good, up to exponentially 



small corrections, so long as the various kinks in each pseudoenergy ea\0) were clearly 
separated along the ^-axis. These kinks have a size of order one (the precise value will 
depend on the model, but in any case we are only interested in orders of magnitude 
in comparison with x and ^0 here), and so the crossover will start when the expected 
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positions of two kinks become closer than this. There are generally two different interkink 
separations for any given (x,^o) (this is rather clear from looking at figure 1); they are: 

Bp - e^^+i)k+s-p = -2a; - ((n-l)/c + 3)60 ; 
9(^ri+i)k+s-p - Op+k = 2x + {nk + s)0o . 

(That these two are positive follows from (|3.1| ); strictly speaking there are also steadily 
growing separations for x larger than the final crossover, after the two spirals have become 
completely disentangled.) The interkink separations are therefore overly small only in 
regions of order one about each crossover value of x. Since the intervals between these 
crossover values grow linearly with 9q, the staircase-like nature of c(x, ^0) soon becomes 
pronounced. This is the evidence for the previously-advertised roaming RG trajectories, 
and from the values of c^, the set of fixed points approached by any particular fiow can 
be read off. The ^('^^ x^«/^('=+^) coset models can be imagined to be located on a grid, 
giving the following skeleton for the large-^o pattern of the k, s (s^O) fiow: 

i 

c{k,2k + s) 
i 

c{k, k + s) 
i 

c{k — s, s) ^ c{k, s) 
The set-up for s—0 is less unexpected: 

i 

c{k,2k) 
i 

c{k, k) 
i 

(massive) 

Before attempting to interpret these results, it is worth seeing that they stand up to 
numerical verification. 



(4.4) 



(4.5) 
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5. Numerical work 



The above has rehed rather heavily on assumptions about the behaviour of the so- 
lutions to ( |2.2|) - the existence of kinks and so on - which, though well-motivated, have 
most certainly not been rigourously derived. Therefore it is worth subjecting the proposal 
to some independent checks, and for this we have solved the equations (|2.2|) numerically 
in a number of cases, discretising the 9 axis at intervals of 0.2, and then iterating ( |2.2| ) 
until c{x, 9o) relaxed to a steady value. To gain five-digit precision, ample for the purposes 
of graph-plotting, typically took from 25 to 30 steps. We have only looked at the case 
Q = Ai - for higher-rank algebras the iteration of TBA equations, even in the usual cases, 
is more tricky - but previous numerical work on higher-rank staircase models P, p!6l , |T7 | 
gives no reason to expect any unpleasant surprises. In particular, in refs. Martins 
proposed and investigated numerically the /c = 2, s = instance of ( p.2| ) with Q — and 
the predictions made above for this case are consistent with his findings. 

Figures 2a and 2b show numerical results for all values of s at /c = 2 and /c = 3, 
respectively. In both cases 6q was fixed at 40, and it can be seen that the agreement with 
the predictions of the last two sections is excellent, even including the final 'corner' of 
( [4.4| ) for s 7^ 0. One amusing feature is that for /c = 3, the two fiows with non-zero s have 
the same infra-red central charge, despite their very different behaviours at intermediate 
scales. At given /c, each pair of fiows s, k—s has this property - at the simplest level just 
a refiection of the fact that c(^('=-^) x^^^V^^''^) = c(^(^) x^('=-^)/^('=)). Note though that 
this is only an equality of central charges - since our system in its current form only traces 
the evolution of the ground-state energy, we cannot distinguish between different modular 
invariants having the same central charge. 



6. Discussion 

The convergence of analytical and numerical results leaves little doubt that the be- 



haviour of the solutions to (|2.2|) is as claimed above, even if this has not been rigourously 
proved. It is then very tempting to suppose that (|2.2|) is indeed the TBA system for 
some relativistic field theory, the exact ground-state scaling function of which is given by 
c(x, ^o)- In the absence of any concrete proposals as to what this model might be, we can 
at least discuss some of its expected properties, and decide whether they are consistent 
from other points of view. 

The first of these properties is that as x increases from the deep ultraviolet, c{x,6q) 
does indeed run through the sequence of numbers Cn given by equation , these being 

a subset of the ^(^^ x^(')/^(^+^) central charges. The function pauses near each for a 
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'time' approximately given by k9o/2 for large Oq, before making a sharp transition to c„_i 
over an interval with size of order one, in the process of which / decreases by k. This is just 
the behaviour needed for a staircase model based on these cosets, being consistent with 
the existence of a family of roaming flows which in the large-^o limit comes to approximate 
the known single-hop trajectories. 

More interesting is the predicted behaviour towards the infrared, as x increases and 
/ becomes smaller. The hopping-by-A; cannot continue indefinitely, since / cannot be neg- 
ative. For s=0, the situation is simple: after a last pause near l=k, the final fiow is to 
c = 0, implying that the theory becomes massive at long distances. All previously-studied 
staircase models have shared this property, which is furthermore in accord with the one- 
hop behaviour predicted for the ^id.id.adj -perturbation of a Q^'^^ xQ^'^^ /Q^'^'^^ coset model. 
To understand the curious final hop that all the staircase systems with s^O undergo, we 
first recall the general small-Z pattern of the interpolating trajectories . Even without 
the TBA, it is clear that once / has become smaller than /c, it is no longer possible for / to 
decrease by a further k. However the perturbation 4>id,id,adj is symmetrical with respect to 
the k and / in the numerator of the coset Q^'^^ xQ^^^ /Q^'^^'-\ and so as soon as / decreases 
below /c, it is natural to expect that the fiow induced by this perturbation now hops by 
decreasing k by an amount /, instead of Z by /c. Such a change in direction for the interpo- 
lating trajectories when / becomes smaller than k is indeed predicted by Zamolodchikov's 
TBA analysis 0. 

For the superconformal discrete series {Q=Ai and /c=2), such a phenomenon was 
discussed in ref. , and can be visualised by referring back to figure 2a, which shows the 
form of the two fiows predicted by (|2.2|) for this case. The series of models with / odd 
continues to hop down in steps of 2 until 1=1, the tricritical Ising model, is reached. All 
of the models up to and including this point have been A^=l supersymmetric, a symmetry 
respected by their 4>id,id,adj perturbations - the field is in fact the top component of a 
supermultiplet, all other components vanishing anyway on integration over the Grassman 
directions. However while the 4>id,id,adj operator of the tricritical Ising model again respects 
the A^=l supersymmetry, it is also the ^13 operator in the sense of the minimal c<l 
series, and as such induces a fiow down to the Ising model, which does not possess such a 
symmetry. The interpretation given in ref. is that the supersymmetry is spontaneously 
broken along this last trajectory. At large 6q we can expect the /c=2, s=l staircase model 
to show similar behaviour, with an approximate supersymmetry at short distances being 
broken near x=6q/2 as the trajectory approximates the final step down to the Ising model. 
Thus a change in the hop direction, from decreasing / to decreasing /c, is not ruled out, 
and should probably be associated with the spontaneous breaking of whatever higher 
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symmetry is associated with a given series. For k=2, this was A^=l supersymmetry, and 
the situation was reasonably under control; for large values of k, or algebras Q other than 
Ai, the situation is less clear, and the discussion much more speculative. The possible 
forms of the symmetries associated with such theories were first discussed in refs. ||2^] , and 
later in ref. [2^ where they were called 'fractional supersymmetries', but in particular the 
lack of a satisfactory generalisation of the Grassman variables to these situations makes life 
difficult. A better understanding of all this will be needed before any detailed implications 
can be drawn for the symmetries of the staircase models. 

To return to the general form of the staircase flows, it is tempting to imagine that 
they should ultimately approximate a complete sequence of these single hops, joined up 
nose-to-tail to form a zig-zag series of flows with the left and right coset indices taking 
turns to decrease, until a diagonal coset {k=l) is reached, at which point the flow would 
be to a massive model. However this is not what happens for the spiral staircase flows 
predicted by ( |2.2D and illustrated in (|4.4|): after only a single hop in a new direction, the 
flow grinds to a halt, the flnal destination in the infrared being the xQ^^^ /Q^'^^ coset. 

In terms of the higher symmetries touched on above, the staircase continues to mimic the 
sequence of single-hop flows only so long as the perturbing operator respects the symmetry 
associated with the k^^ series of cosets - after the final fiow in which the symmetry is 
spontaneously broken, there are no longer any directions respecting this symmetry, and 
the fiow is 'trapped'. 

For an alternative understanding of why the fiow might come to such an abrupt stop, 
first recall some work by Lassig on Zamolodchikov's original staircase model. The 
hopping fiows mimicked in this case were those between minimal models, Aip A4p-i, 
induced by the massless (f)i3 perturbation. Lassig pointed out that at large p an RG fiow 
with properties identical to that predicted by Zamolodchikov can be uncovered within 
a perturbative treatment, simply by adding the slightly irrelevant operator ^31 to the 
original perturbation by the just-relevant operator (pi^. Normally one ignores irrelevant 
operators as having no effect on the infrared destination of an RG trajectory, but this may 
not be correct if the situation under discussion involves a crossover in critical behaviour, 
as here. The subtlety is that even an irrelevant operator with respect to the first fixed 
point may induce couplings to relevant operators near the second, repelling the trajectory 
and sending it on to some new destination of even lower criticality. A detailed analysis of 
the relevant RG equations, to lowest order in 1/p, shows that this is indeed what happens 
for the combined (pi3,(p3i fiows in the minimal series - in fact, both operators are nearly 
marginal and mix under the RG fiow, so that, crudely speaking, the irrelevant (psi operator 
near A4p becomes (pis as M.p-i is approached, and sends the fiow leapfrogging on down 
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the minimal sequence. Thus even an irrelevant operator can cause an interpolating flow to 
change its IR behaviour, a judicious choice managing to replicate itself at the next step, 
and ultimately producing a flow with the characteristic staircase-like form. It is important 
here that the original perturbing operator for the interpolating flow, ^13, flowed precisely 
to (^31 in the infrared - had it flowed to anything else, the subspace of couplings under 
consideration would have had to be enlarged at least to include that operator, spoiling the 
simple leapfrogging picture. To repeat this analysis for a level-A; staircase, it is therefore 
natural to imagine a Q^^^ xQ'^^^ / Q'^^^^^ model perturbed by a linear combination of the 
relevant operator (t)\]y = 4'id,id,adj, which on its own induces the interpolating flow to 
the next model down, and the irrelevant operator (t>\']^^~^'^\ along which the trajectory 
arriving from the model one step higher arrives. (There is also an incoming flow </>/7^''''' 
forming part of a level-/ sequence of hops, but we won't worry about this for now.) In 
ref. 0, Zamolodchikov found the incoming operator in the case Q = Ai to have conformal 
dimension 

A,« = l + ^. (6.1) 

corresponding to the fleld (f>id,adj,id- The flow into 1=1 is exceptional, in that there is no 
operator in the infrared model with this conformal dimension - the adjoint is not among 
the level-one representations of SU{2). Instead, the prediction is that A/fj=2, leading to 
the expectation that here the attracting operator is TT. Staying with this case, which 
for the spiral staircase corresponds to setting s=l, we should examine the effect of adding 
some TT to the (/)i3-perturbation of the (/c, 1) Ai coset model. For large k a treatment 
perturbative in should be valid. Note that even in this limit the operator TT retains 
a scaling dimension of 4, a feature which already distinguishes this case from that of the 
combined ^13, ^31 perturbations. But more important is the fact that TT is a symmetrical 
descendant of the identity, so that the operators generated in repeated operator product 
expansions of ^13 and TT only consist of left-right symmetric descendants of ^13. Since ^13 
itself is already nearly marginal in the large-/c limit, any such descendants will be strongly 
irrelevant. Any operator mixing induced through the RG flow is restricted to operators in 
this subalgebra, and therefore should not affect the flnal destination of the flow. Near to 
the (/c, 1) flxed point there are thus (at least) two different subspaces pertinent 

to staircases: one spanned by ^13 and (/)3i, and one by ^13 and TT. The flrst of these 
is relevant to the original {k=l) staircase of Zamolodchikov, as shown by Lassig, while 
we expect the second to contain the tail-end of the k^^ Ai spiral staircase model at s=l, 
undergoing its flnal crossover before flowing down to the coset model at (/c— 1, 1) in the 
far infrared. The more detailed analysis made by Lassig also gave signs of other flows 
complementary to the staircase, obtained by varying the signs of the couplings. It would 
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be interesting to investigate these possibilities, but for the moment we content ourselves 
with this plausibility argument for the form of the final step of the s=l spiral staircase. 

There do not appear to be any obstructions in principle to applying analogous argu- 
ments in the cases s>l and/or Q^A\. Zamolodchikov found that in the case relevant for 
our Q = Ai, s=2 staircases, the incoming direction at the end of the second-last crossover 
mimicked by the roaming flow is GG, rather than the TT above, where G is the spin 3/2 
part of the supercurrent. This is certainly suggestive, but to give a satisfactory analysis 
of the general case requires a more detailed knowledge of the fractional supersymmetry 
algebras than we have at present. 

A final area for further speculation concerns the field theories underlying our TBA 
systems. If such models exist, they must have a rich and varied structure, and only a 
subset of them can be massive. On the basis of the magnonic structure, we would expect 
to find non-diagonal scattering in these cases. The k=l systems were intimately related to 
the real-coupling affine Toda theories [p|,p!0[| , and so it is possible that further insight may 
be found in the fractional supersymmetric sine-Gordon models introduced by Bernard and 



LeClair [27]. While we have no definite suggestions to make in this direction, it seems to 



be an interesting area for further work. 
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Figure 1 : The arrangement of kinks around tine cylinder. 




Figure 2 : Plots of ground state scaling functions. 
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Fig. 2b 
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